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Abstract 

We extend the phase space slicing method to allow for heavy quarks and 
fragmentation functions. The method can be used to calculate differential 
cross section in which any particular particle (massive or massless) is 
tagged. 
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A next-to-leading order calculation has well known advantages compared to 
a leading order calculation: it reduces the normalization uncertainty, it starts to 
reconstruct the parton shower, and it tests the convergence of the QCD pertur- 
bative expansion. Over the last few years the emphasis has been on constructing 
Monte-Carlo programs that include all a s corrections and are fully differential 
in the final state particle momenta, such that any experimental cuts can be im- 
posed. One such method is the so-called "phase space slicing" method. The for- 
malism for this method developed by Giele, Glover and Kosower Q introduced 
a high degree of automation by using the following ingredients: decomposition 
of the amplitude according to the color structure into colorless subamplitudes, 
factorization of the phase space and of the colorless subamplitudes in the soft 
and collinear region, and the generalization of crossing to NLO. We have ex- 
tended this particular formalism to allow for heavy quarks and fragmentation 
functions. In this short contribution, we briefly review the different aspects of 
the method, full details will be presented elsewhere {1. The method has been 
used in Ref. || to calculate the QCD corrections to W+ heavy quark production 
at the Tevatron. 

The QCD corrections consist of virtual and real corrections. The virtual 
corrections are the interference between the LO and all the one loop diagrams 
and must be calculated in n = 4 + e dimension in order to regularize the singu- 
larities. Coupling constant and mass renormalization take care of the ultraviolet 
singularities, and require the introduction of the renormalization scale. At the 
end, some collinear and soft singularities remain as 1/e and 1/e 2 poles. The real 
corrections are those contributions with one more parton than the LO and have 
soft and/or collinear singularities. One begins by considering the processes were 
all the quarks and gluons are in the final state, e.g., V — ► qq~ + n gluons, where 
V stands for an electroweak gauge boson, such that all the singularities cancel 
without having to do mass factorization. The basic idea of the phase space 
slicing method is to separate the phase space in two regions using the invariants 
Stj = 2Pi.Pj, where the Pi are the momentum of the final state particles. The 
hard region is defined so that all the Sij are bigger than a theoretical cut-off 
S m in- In this case, the calculation can be done numerically in n = 4 dimensions. 
The collinear and soft region is defined such that one or two Sij are smaller than 
S m i n . In this case, the calculation must be done analytically in n = 4 + e di- 
mensions. If S m in is small enough, the soft and collinear approximation can 
be used such that the integration in n dimensions is greatly simplify. We have 
generalized the soft approximation to the case where the particles involved are 
massive. In the collinear region, the mass regularizes the singularities and the 
calculation can be done numerically. The 1/e and 1/e 2 poles that remain after 
the integration over the soft and collinear region cancel with the corresponding 
poles of the virtual contributions. At the end we are left with a "K" factor, 
proportional to the born cross section and dependent on S m in, the finite part of 
the virtual contribution, and the real corrections in the hard region, that also 
depends on S m in- An important numerical test is that any observables should 
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be independent of S m in- 

The generalization of crossing to NLO is done through the use of the so- 
called "crossing functions" . Let us consider the process: pp — > V + n jets. The 
basic idea here is that we do not want to redo the cancellation of all the poles, 
but rather use the "K" factor already derived for V — > (n + 2) jets. First, the 
usual crossing for all the matrix elements is done, along with the crossing of 
the finite terms of the virtual piece and the crossing of the "K" factor with 
appropriate analytical continuation. Then two corrections must be applied in 
n = 4 + e dimensions: 1) subtraction of some collinear singularities included in 
"K" that are not present when the born is crossed, 2) add some initial state 
collinear singularities not yet included in "K" . Along with mass factorization, 
all these corrections give terms that are proportional to the Born cross section: 

a « E/ dxi J dx 2 C^ scheme (x 1 )f^(x 2 )d^ b °(x u x 2 ) . (1) 

a. b 

fa is the distribution function of parton a inside of the hadron H, and the 
are the factorization scheme dependent crossing functions: 

c H,scheme { ^ fa = A H {x) j^fnnnj + ^scheme^ (2) 

where \xp is the factorization scale. The crossing function are universal and 
only need to be calculated once for a given set of parton distribution functions. 

and are convolution integrals of splitting-like functions with different 
parton distribution functions. 

To add fragmentation functions to the formalism we adopted the same idea 
as in the crossing case: we want to use the "K" factor already calculated. Let 
us consider the process: V — > H + (n — 1) jets. First the NLO calculation 
for V — > n jets is convoluted with all the appropriate fragmentation functions. 
Then, corrections must be applied in n = 4 + e dimensions. When the collinear 
contribution from a parton h splitting to parton i and j is calculated for the 
"K" factor, the phase space is not only integrated over Sij up to S m j n , but also 
over the momentum fraction of i compared to h. This last information is needed 
to properly add the fragmentation functions and should not be integrated over. 
Furthermore, this contribution is convoluted with the parton h fragmentation 
function, instead of i or j. The corrections for both of these effects are, along 
with mass factorization, proportional to the Born cross section: 

a s ^ j dzdai T^ scheme {z). (3) 

h 

The are the factorization scheme dependent "tagging functions" , they have 
the same properties and functional S m i n dependence as the crossing functions. 
In the massive case, there is no collinear contribution included in the "K" factor 
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of the calculation without the fragmentation functions. We derived the heavy 
quark tagging and crossing functions that implement the variable flavor num- 
ber scheme Q (it defines the heavy quark factorization scheme). This takes 
care of large logarithms involving the heavy quark mass, leads to collinear safe 
quantities in the sense that when the mass of the heavy quark tends to zero 
the massless result is recovered, and makes our formalism applicable at any 
transverse energy. 
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